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We introduce dissipative effects in the effective field theory of hydrodynamics. We do this in a 
model-independent fashion by coupling the long-distance degrees of freedom explicitly kept in the 
effective field theory to a generic sector that "lives in the fluid" , which corresponds physically to the 
microscopic constituents of the fluid. At linear order in perturbations, the symmetries, the derivative 
expansion, and the assumption that this microscopic sector is thermalized, allow us to characterize 
the leading dissipative effects at low frequencies via three parameters only, which correspond to bulk 
viscosity, shear viscosity, and — in the presence of a conserved charge — heat conduction. Using our 
methods we re-derive the Kubo relations for these transport coefficients. 



I. INTRODUCTION 

Hydrodynamics has recently been recast into an ef- 
fective field theory (EFT) language, with an emphasis 
on its internal and spacetime symmetries, their sponta- 
neous breaking pattern, the associated Goldstone exci- 
tations and their interactions, the derivative expansion, 
and in general on the systematics of the EFT program 
[THS]. For certain questions, this approach offers a num- 
ber of advantages over the traditional one, but so far 
it has neglected a crucial feature of real-world hydrody- 
namics: the presence of dissipative effects. Dissipation 
appears in the gradient expansion of hydrodynamics as a 
first order correction to the perfect fluid equations, which 
are the continuity and the (relativistic generalization of 
the) Euler equations. In fact, for fluids that do not carry 
anomalous charges , all first order corrections are dis- 
sipative. It is then clear that for the EFT program to be 
useful beyond zeroth order in the derivative expansion, 
one has to find a way to accommodate dissipative effects. 

In the standard parameterization, first order dissipa- 
tive effects are characterized by three coefficients: bulk 
viscosity, shear viscosity, and heat conduction. Physi- 
cally, at least for weakly coupled fluids, dissipation arises 
because of the fluid's microscopic constituents' diffusion, 
which tends to erase any gradients the macroscopic quan- 
tities like temperature, velocity field, etc., might have 
. This process effectively converts the mechanical en- 
ergy carried by long wavelength perturbations, like sound 
waves for instance, into thermal energy. Since diffusion 
is essentially unobstructed when the microscopic con- 
stituents are weakly interacting, very weakly coupled flu- 
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1 Bulk viscosity has a different physical origin. See e.g. [8], 



ids are, from the viewpoint of their long-distance hydro- 
dynamical description, the most dissipative ones. This 
is the usual counterintuitive property of shear viscosity: 
at lowest order in perturbation theory, it grows linearly 
with the mean free path (see |9j for a discussion about 
this point.) One could think that going in the opposite 
direction — to very strong coupling — might make diffusion 
and thus dissipation completely unimportant. But, ap- 
parently, this is not case. It has been conjectured that 
there exists an absolute lower bound on the ratio of shear 
viscosity (77) to entropy density (s) [TU1 ITT] . 



Interestingly, heavy ion collision data indicate that the 
quark-gluon plasma has an 77 to s ratio of the same order 
as the proposed bound. 

Systems violating such a bound have been proposed — 
see [12] for a recent review — but it is still an open 
question whether there exists a fundamental bound that 
is just somewhat lower than ([T]). For example, in [T3] it 
has been argued that Eq. ([I]) can be violated for theories 
with gravitational duals, but still there is a (somewhat 
weaker) bound: 2 > which follows after enforcing 

causality in the bulk or micro-causality in the boundary 
CFT. In any case, the mere existence of such a bound 
still defies a purely field-theoretic justification. There- 
fore, part of our motivation to characterize dissipation 
in hydrodynamics in an EFT language, is to derive — if it 
exists — a fundamental bound from sacred properties of 
relativistic quantum field theory. For instance, it might 
follow from unitarity, via dispersion relations [141 115) . 
At this stage we make no progress in this particular 
direction, and keep it open for future investigation. 

Without further ado, we now discuss how to include 
dissipative effects in the EFT formulation of hydrody- 
namics. We will paraphrase a method we learned from 
|16| where absorptive phenomena in black hole physics 
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were dealt with in an EFT fashion. (See also [T?HT§] for 
generalizations, [20] for a somewhat different approach, 
and |21) for an application of the same techniques in a 
different context.) 



II. THE GENERAL IDEA 

Clearly, a field theory with a local action is non- 
dissipative by construction 2 . But so is Nature: In any 
physical system, we call 'dissipation' the transfer of en- 
ergy from the degrees of freedom we are interested in (col- 
lectively denoted by <j), in the following) to others which 
we are not keeping track of (collectively denoted by x)j 
either because we are not concerned about them, or be- 
cause describing them is too complicated or impractical. 
So, the best way to approach dissipation from a field the- 
ory viewpoint — at least conceptually — is to keep in mind 
that these additional degrees of freedom should also ap- 
pear in the action of the system. That is, if we were to 
write the full action for and Xj we would have 

S[4>, X ]=S [4>] + S x \x] + S in% [<t>iX] ■ (2) 

So is the action we would write for 4> alone, if we forgot 
about x- &x g° verns the dynamics of x- Sint couples the 
two sectors, and is responsible for exchanging energy be- 
tween them. If we now compute observables involving our 
<j) only, we can detect 'dissipative' effects — corresponding 
to exciting the x degrees of freedom — which cannot be re- 
produced by using So alone. For instance, the S-matrix 
restricted to the (/>-sector is non-unitary whenever pro- 
ducing x-excitations is energetically allowed. 

In the particular case we are interested in, x stands 
for the degrees of freedom of the microscopic constituents 
making up the fluid. For instance, for a weakly coupled, 
non-relativistic fluid made up of massive point-particles, 
X stands for the positions of these particles. On the other 
hand, c/> stands for the collective degrees of freedom, like 
sound waves for instance, which are those explicitly kept 
by the hydrodynamical description 3 . Notice that at all 
times we are dealing with one and the same fluid, and 
its microscopic constituents. The splitting in Eq. ^ is 
one of the key features of the EFT formalism and the 
(emergent) dynamics in the long- wavelength limit. Hy- 
drodynamics is about the dynamics of <f>. 



2 The formal trick of adding an explicit time-dependence to a La- 
grangian to make the energy not conserved — see e.g. 1221 — might 
work to reproduce the desired dissipative equations of motion, 
but (i) is not systematic, i.e. it is not clear what the rules of 
the game are, and has therefore no predictive power, and more 
importantly (ii) does not correspond to the physical origin of 
dissipation, which is that there are additional degrees of freedom 
that have been ignored. 

3 Strictly speaking, to avoid double counting, one should remove 
from the x's the combinations of the individual particle positions 
that make up the <f>'s. 



To illustrate the general idea, in this section we will 
not commit to the hydrodynamical case, nor will we go 
into many details. Rather, we will keep the discussion as 
general and as schematic as possible. We will only as- 
sume that the interaction Lagrangian S[ n t can be treated 
as a small perturbation. If this is not the case — if the 
two sectors are strongly coupled to each other — then it is 
not even clear how to talk separately of the 0-sector and 
of the x-sector. In other words, we are assuming that as 
a first approximation, one can neglect the x's when talk- 
ing about the </>'s. For hydrodynamics, as will see, this 
will be guaranteed by the symmetries: at low frequen- 
cies and momenta, all the interactions of the 0's become 
negligible, including those with the x' s - Notice that we 
are not assuming anything about interactions within the 
X sector: they can be arbitrarily strong. 

Now, the crucial question is how to make use of ex- 
pression ([2]), without actually specifying what the x's 
and their dynamics really are. The idea is to make the 
dependence of the interaction piece S- m t on </> explicit, 
while keeping that on x implicit. Schematically: 

Si„t = Jd 4 xY / d n m (x)O n , m (x) . (3) 

The O's are 'composite operators' of the x-sector — local 
combinations of the x's and their derivatives. As usual, 
one expects all couplings allowed by symmetry to appear 
in the action. So, in particular, the O's should carry 
spacetime and possibly internal indices in order to make 
the combinations appearing in S lnt invariant under all the 
symmetries that act on the (/>'s. Apart from symmetry, 
as usual in EFT, the other organizational principle in the 
infinite series |3]) is the derivative expansion: terms with 
fewer derivatives acting on the long distance/low energy 
degrees of freedom (</>) matter the most at low energies 
and momenta. 

Now, in any observable that involves measuring the ^'s 
only — like for instance a (</></> •■■(/)} correlation function, 
or a 4>(j> — > 4>4> scattering amplitude — all effects due to the 
presence of the x's, dissipative or otherwise, are "medi- 
ated" by the correlation functions of these O composite 
operators. As an example, consider a coupling (linear in 
(/>) between the two sectors of the form 

S int = \Jd 4 x<f ) 0, (4) 

where A is a small coupling constant. For instance, sup- 
pose that we are interested in computing the T-ordered 
two-point function of <j) in the standard vacuum (i.e. the 
vacuum for both the <p sector and the x sector, for as 
we will see in a moment, computing this same correlator 
with a non-vacuum dissipative x sector will necessarily 
complicate the story). This two-point function will re- 
ceive contributions from Sq and from Sint- We can com- 
pute the latter contribution in perturbation theory for 
A. For instance, if </>'s only interaction is that contained 
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FIG. 1: Feynman diagram representation of eq. the solid 
lines represent the 4> propagators and the gray circles the two- 
point function ofO. 



in iSi nt above, this would correspond to the Feynman di- 
agram series of fig. I. In that case, neglecting combi- 
natoric factors, powers of i, and momentum-conserving 
delta-functions, we would have schematically 



(<I>(p)<K-p)) = (4>(p)<f>(-p)) 

+ A 2 (0(p)0(-p)) o 2 (O(p)O(-p)) 
+ X 4 (^p)H-p))o 3 (0(p)0(-p)) Q 2 



(5) 



where T-ordering is understood, and the subscript zeroes 
denote that those two-point functions are to be computed 
at zeroth order in A, that is, in the absence of any inter- 
actions between (j> and x- Once (4>4>)o an d (OO)q are 
known, the full {4>4>) can be computed at any order in A, 
without any further explicit reference to the x dynamics. 
This is analogous to the standard Feynman-diagram ex- 
pansion for a perturbative QFT, which involves the free 
propagators only. Here the correlators on the r.h.s. are 
not the free ones — they are those determined by So (for 
4>) and by S x (for O) separately. In a more general case, 
where has non-trivial self-interactions and couples to 
the Cs in a more general way, the right-hand side looks 
more complicated because it involves higher-point corre- 
lation functions of <j> and O as well. However, all the cor- 
relators are still evaluated at zero coupling (A) between 
the two sectors. 

As hinted at before, this simple picture gets slightly 
more complicated for correlation functions in more gen- 
eral states and in particular, thermal states. As we will 
discuss at some length in the next section, we will be 
interested in a thermalized x sector. Its real-time cor- 
relation functions and the associated perturbative ex- 
pansion then have to be handled via the so-called In-In, 
or Schwinger-Keldysh, formalism (for extensive reviews, 
see e.g. [23l425j ) . This entails a doubling of the fields 
in the path-integral, 4> — > 4>±, x ~ > X±i which compli- 
cates somewhat the systematics of the Feynman-diagram 
expansion. However, for what we are interested in, we 
can instead consider the effective (linearized) equations 
of motion for the expectation value of <j) that we get by 
"integrating-out" the x sector via In-In path integrals, 
which is essentially an In-In generalization of the quan- 
tum effective action formalism that is appropriate for sys- 
tems described by a density matrix. 

Following the notation of |23j and utilizing the simple 
coupling given by Q , the In-In generating functional for 
the correlation functions of <fi is given schematically by 



J,W\J+,J-} 



const x / T>(j)±Dx± 

c ,i(±S^[0 ± ]±J ± ± ±S x [ X ±]±A0 ± O ± ) 



(6) 
(7) 



where the functional integral over x+ an d X- 1S under- 
stood to include a (thermal) density matrix p(xaiXa) 
for the initial conditions, which are also integrated over 
[55] . As we will see in a second, we will not need to 
be explicit about this. 

Let's assume that (O) = and confine ourselves to 
quadratic order in the <p± fields. Noticing that, from the 
viewpoint of the x sector, the <p± fields act as external 
sources for the operators 0±, we can formally perform 
the functional integration over x+ an d X- an d obtain 



iAl ja n ab ,b 



e iW[J + ,J„] = congtx / V(p±e i(±s 2 [4, ± ]±J ± 4>±) e ^4> a Go' 



where S 2 is the quadratic action for <f>, <j) a = (</>+, 
and Go is a matrix of OO correlators [25] : 



(8) 



( (TO(x 1 )0(x 2 )) (0{x x )0{x 2 )) \ 
yo{xi,x 2 )-^ {0{x2)0{xi)) {0 (x 1 )0(x 2 )T) ) w 

(the T to the right of a sequence of operators implies anti- 
time ordering.) These correlators have to be understood 
as traces involving the density matrix that is appropriate 
for the x sector. 

The In-In effective action T [4>+, <j>—\ is then just the 
Legendre transform of the In-In generating functional 
VF[J + , J~], from which the effective equations of motion 
for ((f)) follow simply as [55] 



ST 



= 



(10) 



However, since we are working at quadratic order in cj>±, 
the effective action T is just whatever appears at the 

after having set J^ 1 to 



exponent in the path integral 
zero: 



r 



S 2 



A 2 ias-rab ±b 

-r<f> Go<P 



(11) 



where two convolutions are understood for the last term. 
We thus get that the linear equation of motion for the 
expectation value of (f> — which, to keep the notation light, 
we also call 4> — is simply 



SS 2 



+ i\ 2 (OO)r *0 = O 



(12) 



where the second term involves precisely the retarded 
two-point function of O: 

{O^O^R^Bih-t^ilOix^Oixi)]) . (13) 

Note that the above conforms to the expectations of the 
usual "linear-response theory" result. What's nice about 
the In-In formalism is that it allows one to generalize 
such a result to all orders in perturbation theory in a 
systematic fashion. 

Keeping these qualifications in mind, and coming back 
to the main message of this section: For generic S- ln t, in 
order to compute observables that involve the 0's only — 
and in particular the time-evolution of (<f>(x)) — we need 
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not be explicit about the dynamics of the x's. We 'only' 
need the n-point correlation functions of the operators 
the </>'s couple to. Of course, knowing all such correla- 
tors is essentially equivalent to having solved the theory 
defined by S x , which, as we stressed, can be arbitrar- 
ily complicated, strongly coupled, or simply unknown. 
Fortunately, in our particular case of hydrodynamics, we 
are interested in such correlators at very low frequencies 
and very long distances only. Moreover, we can assume 
that the \ sector — whatever it is — is in a state of thermal 
cquilbrium. As we will see, this allows us to parameterize 
the leading low-frequency, long-distance behavior of the 
relevant correlators by three coefficients only. 



III. LOW FREQUENCY, LONG DISTANCE 
BEHAVIOR OF CORRELATORS 



Consider the two-point function for a generic opera- 
tor in the \ sector, {0(x, t)0(x', t')). If, in the absence 
of external perturbations — due for instance to our 0's — , 
the x sector is thermalized, then the average (...) has 
to be interpreted as a thermal trace with a density ma- 
trix p oc e~P H , or, in the presence of a conserved charge, 
p oc e~P( H ~^®\ (We will use a quantum mechanical lan- 
guage, but everything we say applies straightforwardly 
to classical statistical systems as well.) Now, we will 
assume that we identified correctly all the degrees of 
freedom that can propagate at long distances and for 
long times-we called them 0-and that we constructed the 
most general EFT for them, encoded by Sq [</>]. We will 
be more explicit in the next section, but for the moment, 
it suffices to say that these 0's correspond to the degrees 
of freedom traditionally associated with hydrodynamics: 
long-wavelength fluctuations in the energy density, in the 
velocity field, in the charge density, etc. Following the 
traditional language, we have 'hydrodynamic modes' — 
i.e. physical variables with non-trivial long-range, late 
time correlators — for each conserved quantity: energy, 
momentum, charge. It is usually believed that thermal 
equilibrium erases all other information that is not asso- 
ciated with conserved charges. In particular, it is usually 
believed that in a thermal system correlators for quan- 
tities that are not densities for conserved charges decay 
rapidly, faster than any power, at very large distances 
and at very late times — roughly speaking, at distances 
and times larger than the mean free path and the mean 
free time, respectively. 

Following this intuition, we will assume that the %- 
sector only features such rapidly decaying correlators. 
As we will see, this does not imply that it does not fea- 
ture gapless excitations. Indeed: if there were no gapless 
X-excitations, it would not be possible for very low fre- 
quency 4> fields to transfer any energy to the \ sector. 
That is: at frequencies lower than the gap, there would 
be no dissipation whatsoever. Now, if an (OO) corre- 
lator decays faster than any power at large space- and 



time-separations, then its Fourier transform 

G(u,k) = J d 3 xdte i ^ t ~^ (O(x,t)O(0)) , (14) 

is differentiable for real u> and k — infinitely many times — 
at uj = k = 0. In particular, it admits a Taylor expansion 
in powers of co and k about the origin. This means that 
at very low frequencies and momenta, we can parame- 
terize our two-point function by just a few numbers — the 
coefficients of the leading terms in such a Taylor expan- 
sion. 

To develop some physical intuition, it is useful to 
rephrase the above statement in terms of the spectral 
density for the operator O. So far we have been cava- 
lier about the ordering of operators inside the two-point 
function. As pointed out in the last section, we will be 
mostly interested in the retarded two-point function, 



G R (x,t) = e(t)([O(x,t),O(0)} 



(15) 



which describes the causal response of the system to 
external disturbances, in the sense that adding a term 
/ d 3 xJO to the Hamiltonian — where J(x, t) is a given 
external source — triggers a response in the expectation 
value of O 



/OO 
dt'd 3 x'G R {x-x', t-t') J(x", t')+0(J 2 ) 
-oo 



(16) 

(we have assumed that the background expectation value 
of O vanishes, i.e (0(x,t)) j =0 = 0). Its Fourier trans- 
form admits the spectral representation 



G R (w,k) 



_ f +c °du>o 

J-oo K 



lu — loq + ie 



p(u> ,k) , (17) 



where p(coo,k) — the spectral density — is a real, non- 
negative function (for positive loq) that quantifies the 
density of states the system has at energy cj and mo- 
mentum k, weighed by the overlap the operator O has 
with them 4 . 

One is often interested in separating the real and imagi- 
nary parts of Fourier-space correlation functions, because 
they contribute to different phenomena. In particular, 
the dissipative effects we are after will be associated with 
the imaginary part of iG R , which, given the distribu- 
tional identity 



— - — = P— — in 6(x) , 
x + ie x 



(19) 



4 The finite-temperature spectral density is given by 

„( W) fc) = | (l - e-"") (Tre-^)" 1 ^e-^» (18) 

n,m 

X (2n) 4 6(u, + E n -E m )8 z {k + p n -p m ) |(n|O(0,0)|m>| 2 

from which the non-negativity (for positive lu) follows immedi- 
ately. 
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is simply the spectral density: 



lm(iG R (oj,k)) = p{uj,k' 



(20) 



Our discussion following (14) thus implies that the rcctly to sect. |IVB 



which encodes the self-interactions of the <fi sector. We 
refer the reader to refs. [U EJ 0] for the details of the 
derivation, and we quote here the main results only. The 
reader already familiar with the notation can skip di- 



spectral density should be infinitely differentiable for real 
uj and k at uj = k = 0, and that it should admit a low- 
frequency, low-momentum Taylor expansion. Moreover, 
standard arguments (see e.g. [26]) imply that the imag- 
inary part of iGr is odd under uj — » — uj (while the real 
part is even), so that in the Taylor expansion of p we 
only have odd powers of ui. The dependence on k is 
constrained by rotational invariance. If O is a scalar op- 
erator, it has to involve 1, \k\ 2 , \k\ A , If O carries 
a vector index i, the fc-dependence of the tensor spec- 
tral density pij will involve the combinations 5 tJ , k l k° , 
| | 2 , . . . ; And so on for higher rank tensors. Given 
these properties, at very low frequencies and momenta, 
the spectral density of a tensor operator that transforms 
irreducibly under rotations can be parameterized by just 
one number — the first coefficient in its Taylor expansion: 



p(uj, k) ~ Auj x <5 • ■ ■ 8 



uj, k -> , 



(21) 



where 8 ■ ■ ■ 8 stands for the combination of Kronecker- 
deltas with the right symmetries 5 . Notice that A has to 
be positive, because p is positive for positive uj. 

We thus see that the absence of long-range, late-time 
correlations in the x sector does not forbid the existence 
of gapless excitations. These can exist, as long as the 
zero- momentum density of states (i) is a regular contin- 
uum in a neighborhood of uj = 0, and (it) goes to zero 
at zero frequency, at least as fast as uj. For instance, a 
(5-function contribution to the spectral density, peaked at 
uj = 0, is not allowed. This would correspond to a gapless 
'single particle' pole in correlators — i.e. to an excitation 
with a power-law propagator at very long distances and 
at very late times. According to our assumptions above, 
this should be included in the <h sector. 



IV. THE ACTUAL COUPLINGS 

We can now be more specific about the structure of 
Sint- To begin with, let us review what goes into So, 



5 For any operator O of given spin s, there is only one possible 
such combination that can appear in the (OO) correlator. The 
reason is that in the tensor product of two spin s representations, 
the singlet (spin-0) representation appears only once: 



(2s + 1) ® (2s + 1) = 1 3 ( 



(4s + 1) 



(22) 



For instance, if Oij is symmetric and traceless, that is, spin 2, 
its two point function at zeroth order in k has to take the form 



(OijO M ) xS ik S jl + s a s ]k 



(23) 



A. Review: the structure of So 

For an ordinary fluid that does not carry conserved 
charges, the infrared degrees of freedom are nothing but 
the individual positions of its volume elements. For our 
purposes, it is convenient to parameterize them via three 
scalar fields (ft 1 (x, t), with 1—1,2, 3, which give the co- 
moving (or 'Lagrangian') coordinates (p 1 of the volume 
element occupying physical (or 'Eulerian') position x at 
time t. Of course the inverse mapping x(cj) 1 ,t) offers 
a completely equivalent — and perhaps more intuitive — 
description, but we find the former more convenient be- 
cause it makes implementing spacetime symmetries im- 
mediate: the cf) 1 behave as scalars under Poincare trans- 
formations. 

There are also internal symmetries acting on the </> 7 's. 
These take a particularly simple form if we choose the 
comoving coordinates in the following way: for an ho- 
mogeneous fluid at rest, in equilibrium at some given 
pressure po, we identify them with the physical ones: 



I Pa 



Then the internal symmetries are 



b 1 ^cb 1 



o z j<p- 

/ ^0 / + C / (0) 



a = const 
O'j e SO(3) 



dct 



1 



(24) 



(25) 
(26) 

(27) 



The first and second lines encode the homogeneity and 
isotropy of the internal space of the fluid. The third 
differentiates a fluid from an isotropic solid: the dynamics 
of a fluid only care about compressional deformations, 
and are insensitive to a reshuffling of the volume elements 
that does not change their local density. 

At lowest order in the derivative expansion, internal 
and spacetime symmetries force the action to take the 
form, 



(28) 



S Q = / d^xFib) , b=Jdetd i 



where F is a generic function, which characterizes — in 
this language — the equation of state of the fluid (we will 
be more explicit about this below). 



The scalar quantity b defined in ( 28 ) , is a (relativistic) 
measure of the compression level of the fluid. Another 
fundamental object is the vector 



(29) 
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It is identically conserved, 



= 



(30) 



as a result of its e-tensor structure, and is related to b 
and the fluid's four-velocity u M via 



b 2 = -J^ 



(31) 



At lowest order in the derivative expansion, J M can be 
identified with the fluid's entropy current s M , and b with 
the entropy density s [J]. 

This formalism can be straightforwardly extended to 
accommodate a conserved charge carried by the fluid, like 
e.g. baryon number jlj [37]. One first introduces a new 
scalar field ip(x,t) that shifts under the U(l) symmetry 
transformation associated with the charge: 



if) ip + a 



const 



(32) 



Then, in order to describe an ordinary charge-carrying 
fluid — as opposed to a superfluid — one promotes such a 
symmetry to a (//-dependent shift symmetry: 



if) -> tp + aty 1 ) . 



(33) 



This guarantees that the Noether current associated with 
the U(l) transformation (32) aligns with the fluid's four 
velocity, 



r 



(34) 



as befits an ordinary fluid at lowest order in the derivative 
expansion. The lowest order action should now read 



S = d 4 xF(b,y), 



(35) 

where b is the same as above, and y is defined as 

y = u " . (36) 

The stress-energy tensor and the charge current follow 
straightforwardly from the action. They are 



T„ v = {F y y - F b b)Bj} d^ 1 d u dp J + (F- F y y) Vt , v 



f = F y u» , 



(37) 
(38) 



where F's subscripts denote differentiation, and B TJ is 
the inverse of the matrix 



b ij = a^^y 



(39) 



Equating T^ v and j M above with the standard forms they 
take for a perfect fluid, and using thermodynamic identi- 
ties, yields the following "dictionary" between field the- 
ory variables and hydrodynamical ones @]: 



P = F y y-F 
p = F-F b b 
s = b, T=-F b 



n = F y 



(40) 
(41) 
(42) 
(43) 



and m m is still given by eq. (31 1. We now see that the 



function F(b,y) is related to the equation of state: for 
instance, from the first, third, and fourth lines we have 



F(s,p) = n(s, p) p - p(s, p) 



(44) 



The field theory language selects entropy density and 
chemical potential as the natural pair of thermodynamic 
variables to work with. 

In the following, we will consider small perturbations 
about an homogeneous equilibrium configuration, which, 
generalizing ( 24 1 to arbitrary pressure and to the case 



with charge, is described in our language by the field 
configuration 



ipo(x) =yot, 



(45) 



where &o an d yo are the fluid's equilibrium entropy 
density and chemical potential. Such a configuration 
spontaneously breaks many of our symmetries: Lorentz 
boosts, completely; Spatial translations and internal (p 1 - 
shifts (eq. (25)), down to the diagonal combination; Spa- 



tial and internal <fi -rotations (eq. (26)), down to the 



diagonal combination; Internal volume-preserving diffs 
(eq. (27)), completely; Time-translations and internal 



■0-shifts (eq. (32)), down to the diagonal combination; 
(//-dependent V'-diffs (eq. (33)), completely. Associated 
with this spontaneous symmetry breaking pattern there 
are Goldstone excitations, which are simply fluctuations 
of our fields about the equilibrium configuration, 



by\x I + n I ), ^)=yo{t + n°) . (46) 



Not all the 7r's feature propagating wave solutions. In- 
deed, after expanding the Lagrangian (35) to quadratic 



order in fluctuations and diagonalizing it, we get [S] 

C ~ \w^ 2 L - C 2(V • 7f L ) 2 ) + i Wo 4 + \F yy yl {^f , 

where ir^ and 7?y are the longitudinal (curl-free) and 
transverse (divergence-free) components of ir 1 , and 7r° 
is a suitable linear combinations of 7r° and 7?l, 



yyVo 



V-7TL 



(48) 



Moreover, wq = (F y yo 
thalpy density (p + p)o, c 
expression involving various derivatives of F 



F b bo) is the equilibrium en- 
2 is a somewhat complicated 



-which cor- 



responds precisely to the standard (dp/dp) s N — and all 



derivatives of F are computed at the equilibrium values 
y and b . 
We 



thus see from (47) 
Goldstones — the longitudina 



that only one of our 
part of ir 1 — has a standard 
quadratic Lagrangian for a gapless field, and wave so- 
lutions propagating at some finite speed c s , with uj = 
c s k. This Goldstone field corresponds to ordinary sound 
waves. The other Goldstones — t?t and tt° — do not have a 
gradient energy. As a result, at the order we are working, 
they have degenerate dispersion laws, uj = 0. 
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For notational convenience, in the following we will not 
differentiate between capital /, J, . . . indices and lower 
case . . . ones. This is consistent because of the orig- 
inal S , 0(3) spacc x SO (3) internal rotation group, only the 
diagonal SO (3) subgroup is preserved by the background 
(45), and under this subgroup the two sets of indices 



transform in the same way. In particular, the tt 1 excita- 
tions transform as a vector field. 



B. Dissipative couplings: Si n t to linear order 

We are finally in a position to write down the couplings 
of our <p's to the x-sector. There is one physical property 
of the x's that we have not yet been explicit about: in 
a sense that we will try to make precise these degrees of 
freedom "live in the fluid" simply because they "make 
up" the fluid — they are supposed to describe all the de- 
grees of freedom of the fluid's microscopic constituents 
that are not explicitly taken into account by the ifi's. 
This requirement alone should fix their transformation 
properties under all the symmetries that act on the </>'s. 

In what follows we will restrict ourselves to the lowest 
order in the derivative expansion and, more importantly, 
to linear order in the tt fluctuations where, as we shall 
see shortly, the coupling to the x sector can be read off 
from basic properties of Goldstone boson interactions. 
In order to generalize our results to higher orders in the 
Goldstone fields, we would need to apply systematically 
the so-called coset construction to our case. The coset 
construction allows one to write the most general inter- 
actions of Goldstone fields among themselves and with 
other degrees of freedom (x, in our case). One only needs 
to specify the symmetry breaking pattern and the trans- 
formation laws of the x fields under the unbroken sym- 
metries, which in our case are suitably redefined transla- 
tions and rotations. Our physical requirement that the 
X's "live in the fluid" should correspond to very specific 
transformation laws under such symmetries. The gen- 
eral coset construction for internal symmetries has been 
worked out in [28 -f3"U] . It was later generalized to spon- 
taneously broken spacetime symmetries in [3T] (see also 
[3"2l [33] for recent applications). Our case presents yet 
another twist, in that certain spacetime symmetries mix 
with internal ones upon spontaneous breaking, in the 
sense that the unbroken symmetries are specific linear 
combinations of internal and spacetime ones. We leave 
addressing such a systematic construction to future work 
and content ourselves with the correct linearized descrip- 
tion, which we discuss next. 

Let us focus first on the case of a fluid without con- 
served charges. Suppose we start from an equilibrium 
configuration in which our Goldstones tt 1 are set zero, 



4>l{x) = bl'\ I . 

Then, let's turn on a small tt 1 perturbation, 
<f> I (x) = b 1 /3 -(x I +n I (x)) , 



(49) 



(50) 



with very mild spatial gradients and time-derivatives. 
Since tt 1 appears as an addition to x l , this is equivalent 
to performing a small spatial translation of the original 
equilibrium field configuration ( 49 ) , weakly modulated in 
space and time: 



</>o ( f ) -> $ (z, t) = <j>o(x + n(x, t)) 



(51) 



We can now be precise about the meaning of "living in 
the fluid" for the x sector: if the comoving coordinates 
(j) 1 are subjected to a weakly modulated spatial transla- 
tion as in eq. (51 1, the x degrees of freedom undergo the 



same spatial translation. But, following standard Nother 
theorem-type logic, under a modulated spatial transla- 
tion with parameter n(x,t), the x action changes by 



s x [x]^s x [ x ] 



(52) 



where T^ 1 is, by definition, the x sector's contribution to 
the Nother current associated with spatial translations, 
that is, the spatial columns of the x sector's stress-energy 
tensor. Therefore we conclude that, at linear order in tt 1 , 



d 4 xd^T^ 



(no charges.) (53) 



Note that d^T^ 1 ^ and so the above interaction is 
non-trivial, since we are not including in T£ % the tt- 
dependent pieces that are required for conservation of 
the total stress-energy tensor. 

A couple of comments about this expression are in or- 
der. First, the coupling above, while invariant under spa- 
tial translations, rotations, and 7r-shifts, does not seem 



to respect the volume-preserving symmetry of eq. (27) 



At linear order this symmetry requires invariancc under 



Tr(t, x) — > jf(t, x) + e(x) , 



. 



(54) 



Since the e parameters are time- independent, we note 
that the 0-component of eq. (53) does respect the sym- 



metry, whereas the spatial parts do not. At the moment 
we have no satisfactory understanding of this issue, but 
we are confident that ( 53 ) describes the correct linearized 



coupling of tt to the x sector, because, as we will see in 
the next section, it correctly reproduces the first-order 
dissipative effects of hydrodynamics. 

Second, the linear coupling of a Goldstone boson to 
the associated current — which we motivated via our "liv- 
ing in the fluid" logic — is likely a very general feature 
of theories with spontaneously broken symmetries 6 . In 
the Appendix we show that the analog of our coupling 



6 This is not — and has no obvious relation with — the usual state- 
ment that the current for a spontaneously broken symmetry 
interpolates the Goldstone particles, in the sense that given a 
singlc-Goldstonc state \p), one has 



(55) 
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holds for a generic theory with a spontaneously broken 
internal U(l) symmetry, and the logic of that example 
suggests that analogous results should apply for more 
generic (internal) symmetry breaking patterns. For spon- 
taneously broken spacetime symmetries there will be ad- 
ditional subtleties, but ignoring them for the moment, we 
are led to postulate that in the case of a fluid with con- 
served charges the leading order interaction Lagrangian 
will read 



Sin 



-Jd 4 x [d^ 1 Tf + 2/ d^ir j£\ (with charges) , 

(56) 

where ir° is the Goldstone excitation of ip, and the asso- 
ciated yo factor comes directly from its definition (46 1: 



j£ is the current associated with shifts of ip, and turn- 
ing on a ir° field corresponds to shifting ip by yo 71 " • As 
we will see, the second term in ( 56 ) will turn out not to 
be the end of the story in this case — hence the question 
mark — but for the moment, notice that, analogously to 
the first term, the second term is non-trivial {d^jtt ^ 0), 
and not invariant under one of the symmetries, 



n (t,x) -t n°(t,x) +a(x) 



(57) 



which is the linearized version of ( 33 1 . We will now check 



that the first term in ( 56 ) reproduces correctly the first- 



order dissipative phenomena associated with bulk and 
shear viscosity — including the celebrated Kubo relations. 
On the other hand, we will see that in order to model 
heat conduction precisely the second term has to be cor- 
rected, both in its overall coefficient and in its structure. 
We leave deriving these corrections from symmetry con- 
siderations to future work. 



V. REDISCOVERING KUBO RELATIONS 

As advertised in sect. [TTJ we can now compute observ- 
ables that involve our Goldstone excitations, and the x 
sector will contribute indirectly to these observables only 
via the correlators of the composite operators that cou- 
ple to our Goldstones. Since the only couplings that we 
have so far are linear in the Goldstones, the observables 
we are able to compute at this point have to do with the 
free propagation of Goldstone excitations. That is, we 
are able to compute the Goldstone attenuation rates. 



A. Fluid without charges 

Consider first a fluid without conserved charges. Its ex- 
citation spectrum — neglecting dissipative effects — is de- 
scribed by the action (47 1, with the w° part omitted. 
We have a longitudinal mode ttl with uj = c s k, and two 
transverse modes ttt with a degenerate dispersion rela- 
tion uj = 0. Consider now one such excitation propa- 
gating in the fluid. Its coupling to the x sector via the 
interaction ( 53 ) will make it slowly decay away, eventu- 



ally transferring all its initial energy to \ excitations. We 
can compute the rate at which this decay process takes 
place at the level of the classical equations of motion 
for the Goldstones. We could also do the computation at 
the level of Feynman-diagram perturbation theory, which 
would be more in line with our field theoretical approach. 
In particular, since the attenuation rates we are after cor- 
respond to imaginary shifts in the excitations' frequen- 
cies, we should compute the x _me diated corrections to 
the poles of the ir 1 propagator. However, as reviewed in 
sect. [TTJ in the in- in formalism each propagator gets re- 
placed by a 2 x 2 matrix of propagators, which, at least 
for our simple computation, complicates unneccessarily 
the systematics of perturbation theory. 

Following sect. H the linearized eom for tt 1 derived 



from the Goldstone quadratic action ( 47 ) , augmented by 
their interaction with the x sector ( 53 1 , is precisely what 
one would naively expect from linear response theory 7 : 



w (uj 2 tt* - cjft'fc* 7T J ) + iC(%(u, k) n j = 



(58) 



where G % f\ is the retarded two-point function of the com- 
bination that couples to tt 1 in ( 53 1 : 



(59) 



and from now on we will use simply (...) to denote the 
Fourier transforms of retarded two-point functions, eval- 
uated at uj and k. Moreover, it will be understood that 
Gji is evaluated in Fourier space, and its uj, k arguments 
will be omitted. 

In the end we are interested in the imaginary parts 
of the eigenfrequencies of the system, which — at leading 
order in perturbation theory — will be related to the imag- 
inary part of iGr. At this point we could parameterize 
the infrared behavior of Imi • (T£ % T£> S ) as described in 
sect. |III| but, before proceeding let us massage this quan- 
tity a little in order to rewrite it in a form that the reader 
familiar with hydrodynamics will recognize. First, notice 
that according to (20), such a quantity is the spectral 



density of a composite operator (T£ l ) in the x sector. 



This interpolation property implies that the full current has 
terms that are linear in the Golstone field, e.g. for relativistic 
theories = fd^it + . . . Here instead we are focusing on the 
terms in the current that depend on other fields — our \'s — but 
not on it, and we are claiming that, in the Lagrangian, the lin- 
ear coupling of 7r to this other sector involves precisely this tt- 
independent part of the current. 



7 We have assumed, as we did in Section that (T^ l ) 7r= o van- 
ishes. In our formalism, the equilibrium expectation value for 
the fluid's full stress energy tensor is given by \37\ , evaluated at 
the equilibrium configuration | |45| l, that is, it is fully captured by 
the (fi's sector action So- 
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We argued that all local operators in the x sector should 
have very well behaved spectral densities near uj = k = 0, 
at least for real uj and k, with a Taylor expansion start- 
ing as const • uj, and continuing with higher powers of uj 
and k. At low energies and momenta, we are interested 
in just that first term, which we can extract formally by 
taking the nested limit 



uj lim 



A lim (lmi.<Tp7> 



(60) 



Given the regularity of our spectral densities in the in- 
frared, we can take the limits in any order. However, 
taking the limits in the order we have written them al- 
lows us to replace T£ l with the total T pl , which includes 
contributions coming from the Goldstone bosons. The 
reason is that at lowest order in the x -71 " interactions and 
in the derivative expansion, the Golstones' contribution 
to any spectral density is a Dirac-delta peaked at on-shell 



values for uj and k. But the limit in (60 1 carefully dodges 
such on-shell values, both for longitudinal (uj = c s k) and 
for transverse (uj = 0) excitations. At the order we are 
working we thus have 



Im(zG^) 



lim i lim Imz • (T^T 

,.,_>.n u ? . \ 



(61) 



Then, using a standard trick — see e.g. [5] — we can use 
conservation of the full stress-energy tensor to set to zero 
terms in the correlator above that have /i or v equal to 
zero. The reason is that, because of conservation, in 
Fourier space we have the operator equation 



rpOa _ _rpk 
UJ 



(62) 



which yields zero if we take k to zero first, like we are 
doing in the limit above 8 . 
We are thus left with 

Im(iG%) ~ ujk k h ■ lim 1 lim Imi ■ (T kl T l3 ) . (63) 

Following sect. we can now split the stress tensor 
operator T lJ into irreducible representations of the (un- 
broken) rotation group, spin and spin 2, 



= T* 3 + T l 2 3 

rpij 1 rij rpkk 

1 — 3° ' 



T 2 J = T lj - \5» T 



kk 



(64) 
(65) 



8 The manipulations we just performed may seem dangerous: in 
fact, in the last section we insisted that is important that the it 1 
does not couple to the full jf", but only to a non-conserved part 
of it, so that the coupling | |53| l is actually non-trivial. There is no 
contradiction however: the divergence — or the fc M — one needs to 
annihilate the full stress-energy tensor doe s not commute with 
our nested limit, so that the r.h.s. in eq. l|61[l is actually non-zero. 



and parameterize the low-energy behavior of the associ- 
ated spectral densities — in the nested limit we are inter- 
ested in — via two free parameters A Q ^ 2 as 

Im i • (T kl T l J ) ~ A uj ■ 8 kl 8 lj (66) 

Imi • {T kl T l 2 J ) ~A 2 uj- (8 kl 8 lj + 8 kl 8 a - l8 kl 8 lj ) . 

We should also mention that the mixed correlator (T T 2 ) 
vanishes at zero momentum, because of rotational invari- 
ance. 



Plugging these parameterizations into eq. ( 63 ) we get 
Im(iG^) ~ uj k 2 [(A + \A 2 )P l l + A 2 P$] , (67) 
where P^ T are the longitudinal and transverse projectors 



p»3 
L 



k l k> 



T 



8 l] - k l W 



(68) 



The reason it's convenient to split this contribution to 
the 7T* eom as a sum of a longitudinal and a transverse 
part, is that the zeroth-order eom has a similar structure: 



\k*hP 7T J -> [(UJ 2 - C 2 k 2 )Pi 



(69) 



Then, putting everything back into eq. ( 58 ) we get im 



mediately the imaginary parts of the (low-momentum) 
eigenfrequencies : 



Auj l 



(A + U 2 ) . 



2w() 



A 



2 k 2 



w 



(70) 
(71) 



These are the attenuation rates for, respectively, the lon- 
gitudinal and transverse modes. We already see two 
important predictions of our field theoretical approach. 
First, the dissipative nature of the coupling (53): these 



imaginary frequency shifts have the right sign to make 
the Goldstone excitations decay in time, since the posi- 
tivity of Aq^ is guaranteed by the positivity properties 
of any spectral density, as reviewed in sect. |III| Second, 
the attenuation rates scale as k 2 at low momenta, which 
agrees with the standard dissipative hydrodynamics re- 
sults. 

But we can go further. Comparing our attenuation 
rates to the standard ones in the literature — see e.g. [S]— 
we find that our parameters Aq^ correspond to bulk and 
shear viscosity, usually denoted by C and rj: 



C = A , V = A 2 



(72) 



Then, our definitions of Aq^ 2 in eq. (66), match precisely 
the famous Kubo relations for bulk and shear viscosity 
[S]. This is the main result of our paper: an independent 
derivation of the Kubo relations via effective field theory 
techniques. 

From now on, we will refer to the nested limit that 
we used above as 'the Kubo limit', and we will denote 
retarded correlators in that limit by 



uj lim — ( lim (■■•)) 



fe-j-0 



(73) 
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B. Fluid with charges 

We can now extend the same analysis to the case 
of a fluid with conserved charges. However, as antic- 



ipated, the interaction Lagrangian ( 56 1 now does not 
work equally well as for the case without charges. In 
particular, it yields the correct Aut ~ ik 2 scaling for the 
attenuation rates, but it does not reproduce the correct 
numerical factors in the Kubo relations for the corre- 
sponding transport coefficients. We blame this on the 



fact that in order to guess the second term in ( 56 1 , in the 



absence of a "living in the fluid" -type argument we ap- 
plied cavalierly the global symmetry lesson of Appendix 
A] directly to our case, which involves spontaneously bro- 
ken space-time symmetries, and is therefore of a slightly 
more subtle nature. We are confident that the coset con- 
struction will shed light on this issue. For the moment, 
let's see whether there exists a minimal generalization of 
( 56 1 that reproduces the correct physics of dissipation in 



this more general case. 

Let's assume that the \ sector still couples to our Gold- 
stones only via the currents T£ v and j£. Then, by (un- 
broken) rotational and shift-invariance, the possible cou- 
plings at lowest order in derivatives are 



rpOO 



(74) 
(75) 



However, given what we learned above by manipulating 
the (T X T X ) correlators in the Kubo limit, we notice that 
correlators involving T x or j x with a zero index will not 
contribute to the imaginary parts that we are interested 
in. For the purposes of our computation, we can thus 
discard the couplings involving those composite opera- 
tors. Furthermore, among the surviving couplings, the 
first is the only one involving the transverse Goldstones 



7T^: the only other possibility is do"K l j x , but recall that 



at lowest order in the derivative expansion tt % t has a de- 
generate dispersion relation uj = 0, which means that, for 
the transverse modes, such a coupling has to be neglected 
at the order we are working. By rotational invariance, 
the transverse modes cannot mix with either ir l L or 7r°, 
and, given our success above in determining their atten- 
uation rate in the absence of conserved charges, we want 
to keep their couplings to the x sector unaltered, i.e., we 
want djTT l T x l to appear with the same coefficient as in 
(53). For the other couplings, we introduce two arbitrary 
coefficients. 

We thus consider the interaction Lagrangian 

S in t - ~ Jd 4 x [drf T x % + B Our j x + C d^ 1 j x ] , (76) 

and determine the values of B and C by computing the 
attenuation rates and matching these to known results. 
We can focus on the ttl-ttq sector only, since for the trans- 
verse modes we have the same couplings as before, and 
the same analysis applies unaltered. On the other hand, 



when we consider the ttl and ttq equations of motion, we 
are sensitive to new (retarded) correlators, 



eom(7r a ) D iG^fnp 



with 



Im(iG°£) ~ B 2 k l k ] Im« • (fj 3 ) K 



~ -BCujk l k j lmi ■ (fj j ) K 
~ k k k l kWImi- {T ki T l i) K 
+ C 2 uj 2 k l knmi- (j l j j ) K ■ 



(77) 

(78) 
(79) 

(80) 



Following the same logic as in the case without con- 
served charges, we have replaced the x-sector's current 
and stress-tensor with the total ones. We have also used 
that the mixed correlator (j l T jfc ) vanishes at zero mo- 
mentum by rotational invariance. 

According to the general discussion of sect. |III| the 
imaginary parts of the (jj)k and (TT)k correlators have 
to scale as uj. As we already saw, in the TT case this 
matches the Kubo relations that determine bulk and 



shear viscosity — see eqs. (66), (72). In the jj case, there 
is an analogous Kubo relation |34[, which relates the heat 
conductivity x to the coefficient of uj 



Imi-{j i f) K = xT{-%-) 2 uj-5 i i 



(81) 



Putting everything together, we can rewrite the imagi- 
nary parts of the iGr entries as 



lm( lG ™)=xT(-^) 2 B 2 ujk 2 



ImliG 



0L\ 



-XT( 



p+p) 



BCuj 2 k 



l^G L R L )^(C+h)^ 2 + xT(-^-YC 



(82) 
(83) 
(84) 



To find the imaginary shifts of the eigenfrequencies, 
it is convenient to re-express the iGr matrix in the tt°- 
ttl basis that diagonalizes the lowest-order quadratic La- 



grangian (47). Given the definition of tt , eq. (48), and 



the structure of our interaction Lagrangian, eq. ( 76 1 , this 
amounts to just replacing 



C -> C - BA \ , 
in the expressions for iGn above: 

^S6\ _ ,/W n \ 2 r>2 , ,7,2 



~ FyyVO 



Im(zG£) 
Im(zG^) 



,-xT(^) 2 BCuj 2 k + xTB 2 Ak 3 



(85) 

(86) 
(87) 



Im(zG^) = (C + lv) ujk 2 + xT(-^Y (C - BA^ 



We thus get that the eigenfrequencies of the system — 
identified by the vanishing of the determinant of the ma- 
trix defining the eom — get shifted by 



yyVo 



A -6---Xr(-f p ) B 2 



Auj l ~ —i ■ (C+lv)^+xT(^y ( -^-^ 



w cj 
2( 



(89) 



(90) 
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where we have kept only the leading order in k. 

We are now in a position to match our computations 
to the classic dissipative fluid results, which, in the case 
of a fluid with conserved charges, are quite messy — see 
e.g. |35j . At low momenta, the attenuation rates for the 
scalar modes are 9 



Auj. 



(p+p)(dp/dT) r , 



M.und — 1 ' 2(p+p) 
x ((P + P) - 



(C+U) + x(dp/dT)-i 



2T 



dp I 
dT I 



n 2 T dp I 
u s ± dn\T 



n dp I 
cl dn\T 



(91) 
(92) 
A: 2 



After a messy computation (see Appendix |B| one can 
see that, despite their complexity, these expressions agree 
with ours above if we simply choose 



B=-C=- 



yo wo 

b F b 



(93) 



That is, the two coefficients in (761 have in fact the 



same value (up to a sign), which in hydrodynami- 
cal/thermodynamical terms is simply the combination 
fi(p + p) /Ts. The emergence of such a simple final result 
from a long series of fairly messy intermediate steps, gives 
us confidence in the correctness of ( 76 1 with this partic- 



ular choice of coefficients. This completes our matching 
computation. 



VI. DISCUSSION & OUTLOOK 

It is useful to summarize the salient features of our 
results: From purely symmetry arguments and the prin- 
ciples of EFT, we were able to derive that the coupling 
of hydrodynamical modes to a generic thermalized sec- 
tor that "lives in fluid" yields dissipation, with atten- 
uation rates scaling as k 2 . This matches well known 
features of dissipative effects in hydrodynamics, and is 
quite independent of the precise structure of the cou- 
plings that we would write down, following essentially 
from the thermal nature of this extra sector. For flu- 
ids without conserved charges, the living-in-the-fluid re- 
quirement is strong enough to determine — via symmetry 
considerations — the precise structure of the interactions, 
thus allowing us to re-derive Kubo relations. For flu- 
ids with conserved charges, we adopted a "symmetry- 
inspired" ansatz for the interactions, with two free pa- 
rameters, which we determined via matching a proce- 
dure. The emergence of a remarkably simple interaction 
Lagrangian, 



d i x[d^T^ + ^^{d l ^^d^)j' x } , (94) 



9 In fact, ref. |35| computes only the attenuation rate for the sound 
mode. However, following that paper's derivation it is easy to 
spot another scalar mode, which corresponds to our fro and which 
we call the "heat mode" , with an attenuation rate as given below. 



from a long series of fairly cumbersome intermediate 
steps, gives us confidence in the validity of the arguments 
behind our ansatz. 

We feel that these are significant accomplishments. 
However, for our techniques to be a useful tool rather 
than simply an alternative derivation of well-know fea- 
tures of hydrodynamics, we need to extend our analysis 
beyond linear order, and, more importantly, we need to 
understand the systematics of the symmetry structure of 
the dissipative couplings. As we already emphasized, we 
believe that the coset construction will help us in these 
directions. Until then, we are left with a puzzle: Our cou- 
plings to the x sector do not preserve all the symmetries 
we started with, in particular the volume preserving diffs 
(54 1, and the "modulated" shift (57). Because dissipa- 



tion involves time derivatives — the factor of w determin- 
ing the infrared behavior of spectral densities — , and the 
symmetries above are time-independent, this was not an 
issue on the way to reproduce known results. However, 
symmetry breaking terms can be generated, for exam- 
ple if we compute the effects associated with the real 
parts of our correlators, e.g. Re i{T£>T* 1 ), If these have 
frequency-independent pieces, they can yield symmetry- 
violating terms in the Goldstone effective action, like for 
instance a gradient energy for the transverse modes, in 
contradiction with standard properties of hydrodynam- 
ics. Unfortunately, from standard analytic properties of 
retarded Green's functions one can derive a dispersion 
relation of the form (see e.g. [26] ) 



Re(iG fl ( Wo ,0)) 



lm(iGn(uj, 0)) dui 



l-oo w - w 



(95) 



The RHS is strictly negative-definite when we take uj - 
— because of the strict positivity of spectral densities- 
which leads us to the result 



lim Re(iG R (uj)) < 



(96) 



That is, the unwanted frequency-independent pieces are 
in fact forced to be there. One could in principle add 
symmetry-violating local counter-terms to cancel out 
the undesired, symmetry-violating contributions from 
the real parts of our Green's functions (these are ana- 
lytic in frequency and momentum, and therefore local 
in position space). But this would correspond to fine- 
tuning certain Lagrangian terms to zero, which would 
go against the whole point of insisting on symmetries 
as the guiding principle to construct effective field the- 
ories: the only robust properties of physical systems 
should be those ensured by symmetries. A possible res- 
olution to this puzzle may be connected to the validity 



of Eq. (95), which requires for convergence of the RHS, 
lim w _>.o Im (iGji(uj, 0)) — > 0. Even though at this stage 
we are not in a position to argue against such a behav- 
ior, it is nonetheless not required by any basic principle. 

We conclude with an intriguing application of our re- 
sults to black-hole physics. Ref. [H] — from which we 
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borrowed the techniques of sect. [Il] — analyzed dissipative 
effects in the dynamics of black holes, starting from a 
matching computation involving the absorption of gravi- 
tational waves by a black hole in isolation. A classical GR 
computation for such a process yields — in the language 
of sects. [Tl] and |III| — a spectral density for the relevant 
composite operator scaling as oj at low frequencies |16j . 
Our arguments of sect. |III| show that such a behavior 
is characteristic of spectral densities for local operators 
in thermal systems. That is, a black hole absorbs like 
a thermal system. This is yet another indication of the 
thermal nature of black holes. What is remarkable is that 
such an implication here follows from a purely classical 
computation in GR. 
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broken by the vev of <f>: (<f>) 
we can parameterize (f) as 



Then, as it is standard, 



4> = {v + P y*, 



(A4) 



with 7r being the Goldstone boson associated with the 
symmetry breaking, and p the (generically) heavy radial 
excitation, which can be ignored at very low energies. For 
the sake of argument, let us thus set p to zero from now 
on. Note that the symmetry is now realized non-linearly 
on 7r, i.e. 7r — > 7r + a. 

We can expand the action as 



(A5) 



the same way we expanded tfi 1 = x + it 1 for the fluid. 
However, in this parameterization of the fields it is not 
obvious that 7r is derivatively coupled, as guaranteed from 
standard soft-pion theorems for the emission of a single 
soft 7r quantum. For instance, from (A3) we get a cou- 
pling 



C x D \v 2 (e l2q * X * 2 +h-c 



(A6) 



which does not involve derivatives of n. This stems from 
a suboptimal choice of the field variables, and is easily 
remedied via a non-linear redefinition of the x field, which 
as usual does not change the S-matrix: 



X = X e 



■iqn 



(A7) 



Notice that the new %' Re\d is invariant under the [7(1) 
symmetry — the transformation of x is now carried solely 
by the e lq7V factor — and the action becomes: 



C = C 4> [ve™]+C x [x'e 



(A8) 



Appendix A: Linear couplings of a (7(1) Goldstone 

Consider a pair of complex scalar fields charged under 
a [7(1) global symmetry. For consistency of notation with 
our fluid case, let us denote them as <f> and X- Then we 
have a Lagrangian of the sort 



Let's focus on the C x part. By the U (1) symmetry — 
which now only acts on 7r — this action must be in- 
variant under constant 7r shifts, 7r — > 7r + a. Then, 



interpreting the ir in C x [x' e 



T ) as a weakly 



spacetime-dependent [7(1) transformation parameter, 
from Nbther's theorem we get 



c = cm + c x [ x ,4>] , 

invariant under the [7(1) transformation 



X 



e lqa X 



(Al) 



(A2) 



(we are allowing for different charges for <p and X-) No- 
tice that for the moment we are using a slightly different 
notation from the main text: we are including in C x both 
the x-sector's dynamics, and its interactions with the 4> 
sector. For instance, C x might contain interactions of the 
form 



£ X D\( 



(A3) 



As in the case of spatial translations for our fluid, let us 
imagine now that this (7(1) symmetry is spontaneously 



= L 



xlX 



(A9) 



where J£ is the x _sec tor's contribution to the [7(1) 
Nother current, and we omitted terms with more 7r's or 
more derivatives. In other words, at linear order in 7r and 
at lowest order in the derivative expansion, the interac- 
tion between 7r and the x sector has to take the form 



(A10) 



which is exactly the [7(1) analog of our eq. (53). Notice 
that, at this order, it does not matter whether we evalu- 
ate J£ at x or x', since their difference is of first order in 
7r. Notice also that we never really used that x is scalar. 
Clearly our proof is completely general and holds for any 
set of charged fields Xi °f an Y spin. 
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As an alternative, quicker derivation of the same re- 



sult, we can go back to eq. (A5) and use the follow- 



ing trick (a version of Stiickelberg's trick): We promote 
the global symmetry to a gauge symmetry by intro- 
ducing an auxiliary gauge field A^ which transforms as 
Afj. — > A^ — d^a(x), and replace standard derivatives by 
covariant ones. Expanding the action to linear order in 
Afj, we have 



C[<l>,X,A ll ]=C[(l>,x] 



J^ + 0(A 2 ), (All) 



where J^{4>, x] is the conserved current (in the absence 
of _4 M ) for the U(l) global charge. 10 Since we promoted 
this symmetry to a gauge transformation we are guaran- 
teed that 7r disappears from the action, because it can be 
absorbed into A^ by choosing a — — n. This means that, 
to linear order in tt, we wind up with the coupling 11 



(A12) 



where J£ is the x~dependent component of the full 
current at zeroth order in 7r. 12 Hence we conclude that, 
at leading order in the perturbations, the Goldstone 
boson couples to the 7r-independent part of the current 
associated with the broken symmetry. 

The introduction of A^ is equivalent to working in the 
so called unitary gauge, where the Goldstones are set 
to zero and their interactions are encoded in the gauge 
field. The previous analysis suggests that one could per- 
form similar manipulations in the case of our fluid, where 
the Goldstone fields it 1 are associated with the break- 
ing of spatial translations. Now, to go to the unitary 
gauge we must introduce the gauge field associated with 
spatial translations, namely the metric perturbation h^j, 
and the broken generators are the T^ 1 components of 
the stress energy tensor. Hence, the coupling must read: 
T^hni. To introduce the pions we do as before, which 
in our case entails schematically 



This viewpoint might prove useful in extending our re- 
sults to non-linear order. 



Appendix B: Matching the attenuation rates 

In order to match our computed attenuation rates for 
the sound and heat modes given by (|89| and (90) with 
the known values given by (|9 1 h and 



(1921) and fix our free 
parameters B and C we must express many thermody- 
namic quantities such as c 2 , (dp/dT) n , etc. in terms of 
our effective theoretical variables b, y and various deriva- 
tives of F. Explicitly they are as follows: 



dp 
df 

dp 
df 
dp 
dn 



d{F - F b b ) 



d(F y y Q - F) 



(Fy - F b yb ) 2 ~ blFyyF bb 



d(F H /6)=0 



w a F, 



.'/.'/ 



(Bl) 



d(-F + F y y ) 



d(-F b ) 
d(F - F b b ) 



dF v 



V=0 ( F by - F bbFyy) ' 



(B2) 



d(-F b ) 
d{F - F b b ) 



= bn 



FqiF'. 



yfyb 



d{Fy) 



dF v =0 



dF b =0 



FbbF _ F 2 » ( B3 ) 
r bb r yy r by 



F v F bb 



Fhn - FbbFy 



(B4) 



(For c 2 : the F y /b = n/s = const constraint is equivalent 
to the usual one, S — const, N — const [6].) Along 



with the relations (40 )-(43 1 this is all we need to perform 



the matching. We find that for the two computations to 
coincide 



ij = -c wiih m = _yjw = »(p+p) (B5) 



b F b Ts 



h^i = d t 



(A13) as quoted in the text. 
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